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Neutrinoless double beta decay has been the subject of intensive theoretical work as it 
represents the only practical approach to discovering whether neutrinos are Majorana 
particles or not, and whether lepton number is a conserved quantum number. Available 
calculations of matrix elements and phase-space factors are reviewed from the perspective 
of a future large-scale experimental search for Oufif} decay. Somewhat unexpectedly, a 
uniform inverse correlation between phase space and the square of the nuclear matrix 
element emerges. As a consequence, no isotope is either favored or disfavored; all have 
qualitatively the same decay rate per unit mass for any given value of the Majorana 
mass. 
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1. Introduction 

The matter asymmetry of the universe remains one of the deepest mysteries in 
physics. The absence of significant amounts of antimatter requires, as Sakharov 
explained,'^ a time when the universe was not in equihbrium, the non- conservation 
of baryon number, and violation of CP invariance. Non-conservation of baryon and 
lepton number has not been experimentally discovered, despite heroic efforts to ob- 
serve proton decay. The observation of neutrinoless double beta decay would demon- 
strate the non-conservation of lepton number, and by inference the non-conservation 
of baryon number. It is the only known practical approach to discovering whether 
or not neutrinos are their own antiparticles, i.e. 'Majorana' particles. 

Neutrinoless double beta decay can be induced by the exchange of a massive Ma- 
jorana neutrino with an electron-flavor admixture. From neutrino-oscillation data 
neutrinos are known to have mass, and are expected on general theoretical grounds 
to have a Majorana character. The "see-saw" mechanism proposed to explain the 
lightness of neutrinos (see Ref. [H and references contained therein) predicts that 
they are Majorana particles. The effective Majorana mass may be complicated by 
other kinds of new physics and might not be the physical mass of a particle, but it 
serves as a metric for designing and comparing experiments. 
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The rate of neutrinoless double beta decay (Oi/fip) depends on the available 
phase space, the size of the nuclear matrix element, and the effective Majorana 
neutrino niassEE Specifically, the half-life given by 

2 

(1) 



^1/2 



where G^^J is the phase space factor as defined and tabulated by Kotila and 
lachellcP, Mq^, is the nuclear matrix element, and the Majorana mass is: 

(rriee) = I f/|imi + C/^'smse" + U^^mse'^ \ . (2) 

The convention for the Majorana phases is the one given by RodejohannPThe Ma- 
jorana mass is a coherent sum over mass eigenstates with (potentially) CP-violating 
phases, and cancellations can occur. The effective mass can also be modified by in- 
terference with other hypothesized non-standard-model processes. 

Historically the effective axial-vector coupling constant gA has usually been in- 
corporated in the phase space factor, and sometimes in the nuclear matrix element. 
Here, following Kotila and lachcUo, we break it out explicitly. The phase space fac- 
tor g'^J has recently been reevaluated by Kotila and lachello (Table III in Ref. [51) 
with an exact treatment of screening, resulting in significant downward corrections 
of as much as a factor of 2 for the heaviest nuclei. These authors use for the nu- 
clear radius R — vqA^^^ with tq = 1.2 fm, cautioning that some authors have used 
To = 1.1 fm. 

The presence of g\ in the rate introduces a significant uncertainty in the calcu- 
lated rates, in addition to the well-known uncertainty in the nuclear matrix element. 
Barea et al^ and EjirP have fitted the known half-lives for 2z^/3/3 decay and find 
effective values of gA of about 0.8 for shell-model calculations and 0.6 for the Inter- 
acting Boson Model (IBM). Barea et al^also note a weak A-dependence, which we 
neglect here. With these renormalized values there is qualitatively good agreement 
with the data on 2i'/3f3 half-lives. In contrast, the calculated phase-space factors for 
neutrinoless decay are generally presented with the free-nucleon value gA = 1.269, 
gA — 1.25, or gA = 1- The difference between the free-nucleon value for gA and 0.6 
corresponds to a factor of 20 in rate. The extent of the renormalization of gA in 
neutrinoless double beta decay remains a topic of discussion among theorists. 



2. Application of Theory to Experiments 

Experimental work on at least nine different double beta unstable nuclides is in 
progress around the world, not counting the more unconventional projects involving 
electron capture and positron emission. Experiments using large amounts of ^^Ge, 
-'^^^Te, ^"^^Xe, ^^°Nd are actively being prepared or planned by international teams 
involving US researchers. The Super-NEMO Collaboration is focusing its attention 
on a large ^^Se experiment.!^ A fundamental requirement is that, independent of 
technical issues concerning background, resolution, etc., there must be sufficient 
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signal to detect. The sensitivity to a particular value of the Majorana mass depends 
on the phase space for the decay, the nuclear matrix element, the effective value in 
the nuclear medium for the axial- vector coupling constant qa, the equivalent mass 
of isotope M (containing TV atoms) that is actively monitored at 100% efficiency, 
and the time over which the measurement is made. The decay rate per unit mass 
of isotope is 

N hl(2)A^ (0) 4 2|, ,|2 

M Am^ 

= Ho.9\\Mou\^ \{m,,)\^ (3) 

where A is the atomic mass of the isotope, and Na is Avogadro's number. Constants 
are aggregated with g'^J to form the 'specific phase space' Hq^. The phase space 
G^^J is an activity per atom, whereas the specific phase space Hq^, is an activity 
per unit mass. The nuclear matrix elements have been calculated by a number of 
methods, and we make use of a recent compilation by Dueck et alf^ supplemented 
with results for ^^^Te in Refs. [TTl [T^ (for a still more recent evaluation, see Ref. 
[T^ . The matrix elements compiled by Dueck et al. and shown in Table [T] were 



Table 1. Phase-space factors GqJ in units of xlO~^^ from ref. [S] specific phase space Hoi, in units 
of Mg"'^ eV~2 from Eq. [s] and nuclear matrix elements from refs. 1101 1141 for neutrinoless double beta 
decay candidate isotopes. 



Isotope 


Q 




Ho. 


Shell 


GCM 


QRPA 


QRPA 


IBM 


IBM 


PHFB 


PHFB 




keV 




Model 




low 


high 


low 


high 


low 


high 


48Ca 


4272 


24.81 


826.2 


0.85 


2.37 








2.00 






76 Ge 


2039 


2.36 


49.6 


2.81 


4.60 


4.20 


7.24 


4.64 


5.47 






82 Se 


2995 


10.16 


198.1 


2.64 


4.22 


2.94 


6.46 


3.81 


4.41 






96 Zr 


3350 


20.58 


342.7 




5.65 


1.56 


3.12 




2.53 


2.24 


3.46 


lOOMo 


3034 


15.92 


254.5 




5.08 


3.10 


6.07 


3.73 


4.22 


4.71 


7.77 


llOpd 


2018 


4.82 


70.0 












3.62 


5.33 


8.91 


116Cd 


2814 


16.70 


230.1 




4.72 


2.51 


4.52 




2.78 






l24Sn 


2287 


9.04 


116.5 


2.62 


4.81 








3.53 






l28Te 


866 


0.59 


7.4 




4.11 


3.50 


6.16 




4.52 






130rpg 


2527 


14.22 


174.8 


2.65 


5.13 


3.19 


5.50 


3.37 


4.06 


2.99 


5.12 


l36Xe 


2458 


14.58 


171.4 


2.19 


4.20 


1.71 


3.53 




3.35 






i«Nd 


1929 


10.10 


109.1 












1.98 






isONd 


3371 


63.03 


671.7 




1.71 


3.45 




2.32 


2.89 


1.98 


3.70 


I54sm 


1215 


3.02 


31.3 












2.51 






160Gd 


1730 


9.56 


95.5 












3.63 






198pt 


1047 


7.56 


61.0 












1.88 







renormalized by them to a common value qa — 1.25. We renormalize again to 
remove qa from the matrix element entirely, i.e. divide the tabulated numbers by 
1.25^ and square to obtain the quantity |Mo,y|^ of Eq. [i] 

The possible range for the effective value of qa in each case is taken to be not 
larger than the free nucleon value, 1.269, and not smaller than 0.8 for shell-model 
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matrix elements or 0.6 for other calculation methods (GCM, QRPA, IBM, and 
PHFB). The lower 'limits' are guided by the results reported by Barea et aZjZIHand 
based on the renormalization needed to fit experimental lifetimes for 2i//3/3 decays. 
However, there is no consensus that such a large renormalization will apply to Ovf3f3 
decays. It can be argued that in Ov/S/S a closure approximation is applicable, and 
all intermediate states and multipoles are included, whereas in 2j//3/3 the strength is 
concentrated in just a few intermediate l"*" states. Quenching of gA is a symptom of 
basis truncation, and it might therefore be less of an issue for Oi^/3/3. Furthermore, 
many calculations make empirical adjustments to internal parameters such as gpp, 
the particle-particle strength in QRPA, or occupation numbers determined from 
two-nucleon transfer reactions. A detailed discussion is given by Faessler et al^^ 
These adjustments have the effect of correcting for some of the basis truncation, a 
correction that may also diminish the need to renormalize gA- 

While the coupling constant g\ is formally a factor in the nuclear matrix element, 
it has been given an independent role in theoretical calculations. The question of 
how gA is renormalized is decoupled from other sources of uncertainty, such as 
the initial- and final-state wave functions and short-range correlations. To keep 
these two classes of uncertainty distinct, in the following two plots we multiply the 
specific phase space by g\ and plot the product on the vertical axis against the 
square of the nuclear matrix element on the horizontal axis. Since the specific phase 
space itself is essentially exact, this provides a convenient way of displaying both 
sources of uncertainty, and correlations between them can be included on a theory- 
by-theory basis as mentioned. Figure [T] summarizes the rates and theoretical ranges 
for four important isotopes. It may be seen that the four isotopes have comparable 
sensitivity. The vertical span of the regions results from the spread in gA and in a 
Bayesian sense may be thought of as an uncertainty. The true value is likely to be 
found in this range. For the horizontal span, an interpretation as an uncertainty is 
much less satisfactory. Different theoretical approaches do explore different types 
of deficiency in theory, there being no exact theory for this process, but they do 
not necessarily include the true value. Moreover, it can be expected that better, 
or at least different, theories applied in the future might expand the ranges shown: 
no mechanism can reduce a range that encompasses all values. Such behavior is in 
contrast to what one expects for an uncertainty. Despite these shortcomings, there 
is comparative information available at a glance about the theoretical situation for 
each isotope in these plots. 

For clarity we present the results in two plots, the second plot, Fig.[2j containing 
the same isotopes as Fig. [l]but in addition the isotopes ^^Ca, ^^Se, ^^Zr, iooMo, 
-'^^^Pd, ^^®Cd, and ^^^Sn, confining our attention to Q-values > 2 MeV. In this fig- 
ure a correlation between the phase-space factor and the nuclear matrix element is 
strikingly apparent, and there do not even seem to be significant exceptions. This 
seems quite surprising, and to our knowledge has not previously been remarked on. 
The trends in these plots are linear in a log-log display but would be hyperbolae 
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0.1 1 10 100 

Matrix element squared 

Fig. 1. Regions in the renormalized specific phase space g^Hg^ = gjj In(2) ^o'^'^ ^.nd matrix 

element squared |Mo^|^ that encompass modern theoretical calculations, for the candidate neu- 
trinoless double beta decay isotopes '^^Ge, ^^"Te, ^^^Xe, and ^^"Nd. The vertical span refiects 
the range of g^, which differs for the shell model and other models, leading to non-rectangular 
boundaries. The matrix-element calculational methods are shell model (SM), generator-coordinate 
method (GCM), quasiparticle random-phase approximation (QRPA), interacting boson model 
(IBM), and Projected Hartree-Fock Bogoliubov method (PHFB), as given in Table ^ The lines 
indicate the effective Majorana mass that would correspond to a count rate of 1 event per tonne 
per year. 



and not easily seen in linear-scale plots. A different view with the same information 
is shown in Fig. |3] Here the geometric means of the squared matrix elements are 
plotted for a phase-space factor with gA = I- The nucleus ''^Ca is generally con- 
sidered to have a 'hindered' Ov/3/3 matrix element, but it does not appear to be in 
any sense unusual in these figures. One general conclusion is that (at the level of a 
factor of 3 or so), there are no especially favored or disfavored isotopes for a OvPfS 
search. They all have roughly equivalent sensitivity from the theoretical standpoint. 
The specific activity may be expressed as 

N A 2 

Aoi.;^ = OO^SA l(™ee)| (4) 

log (^^oujj^ = log «0i. + 4 log -I- 2 log I (mee) I , 

where gq^ = {Hq,^\Mq^\'^) is a constant with a value ^ X0^'^*° '^ decays per year per 
tonne per eV^. Of course, for experiments that are potentially very costly, factors 
of a few are important, and the deviations from this general trend require the best 
theoretical treatment possible. 

Why would a large phase-space factor imply a small matrix element and con- 
versely? For any given choice of gA, the specific phase space as defined above depends 
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2 3 4 S 6 7 8 9' Z 3 4 5 6 7 8 9' 2 3 4 5 6 7 8 9 

0.1 1 10 100 

Matrix element squared 



Fig. 2. As Fig.[T]but with the addition of the isotopes "SCa (2.2, 2143), ^^Se (17, 514), ^^Zr (13, 
889), if'OMo (25, 660), ""Pd (33,181), "''Cd (9, 597), and ^^*Sn (10, 302). The number pairs are 
the coordinates of the upper rightmost corner of each area, in lieu of labehng. It is more difficult 
to see the details but the overall trend of a correlation between the phase space factor and the 
square of the nuclear matrix element is brought out. 




0.1 1 10 100 

Matrix element squared 



Fig. 3. For each candidate isotope a point is plotted at the geometric mean of the squared matrix 
element range limits (as shown in Fig. [2jl and the phase-space factor evaluated at gA=^- The 
points in order of increasing abscissa value are: ''**Ca, ^^'-'''Nd, ^^^Xe, ^^Zr, "^^Cd, ^^'^Sn, ^^°Te, 
82Se, ■^^Ge, iooMo, and "OPd. 



explicitly on A, and implicitly on the Q-value and Z. The presence of the factor A 
in the specific phase space has little effect on the correlation but slightly reduces 
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Fig. 4. The matrix elements determined individually theory by theory and listed in Tableware 
plotted against the specific phase space. Where a theory predicts a range of values, the geometric 
means of the high and low predictions are plotted. Red triangles: shell model, green circles: GCM, 
blue squares: QRPA, black diamonds: IBM, red crosses: PHFB. 



the scatter. Most of the correlation therefore must be traceable to the Q-value and 
Z. A natural question is whether individual theoretical predictions show the same 
trends as do the summary data. The points are shown in Fig. [4] and include isotopes 
with Q < 2 MeV. The trend is less apparent, with the scatter of values from each 
theory dominating the distribution and the low-Q cases falling below the rest. That 
is to be expected because a power-law scaling relation can only apply in a limited 
range of Q since, while the phase space factor can be arbitrarily small, the matrix 
element cannot increase without limit. 

In Fig. [5] a histogram of the individual entries for the relationship given in Eq. [4] 
is shown, for isotopes with Q > 2 MeV, and with = 1- The distribution shown 
in this figure is the basis for the numerical constant in Eq. |4j Even though the 
individual priors for the matrix elements are not known, the central limit theorem 
evidently leads to a distribution that is approximately Gaussian. Consistency with 
the central limit theorem over all isotopes and all theories lends support to the con- 
jecture that the various theoretical approaches are all calculating a single common 
observable aoi/, within the uncertainty expressed in the width, and the best estimate 
for that observable is the mean of the normal distribution. 

A reasonable estimate for the range of possible values for is 0.8 < gA ^ 1.269. 
The lower limit of the range is 0.8 instead of 0.6 in recognition of the correlation 
between the theory used to calculate the matrix elements and the extent of renor- 
malization of required. (This effect is responsible for the shape of the lower 
boundaries of the areas shown in Figs. l]2 ) If the prior for is assumed to be flat 
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Fig. 5. The individual values of the specific activity determined by theory and listed in Table ^ 
are histogrammed for qa = 1- No averaging is applied in this histogram; for theories with high 
and low entries, both are included. 



in log the effect on the width ol the distribution is relatively modest. Including 
this contribution in quadrature increases the standard deviation from 0.5 to 0.6: 



ao^ = 10^ 



In the future, if a more detailed understanding of the appropriate renormalization 
of gA for each theory emerges, there is the possibility that the width of the distribu- 
tion could actually be reduced, perhaps quite substantially. Simply using = 0.8 
for the shell-model entries and 0.6 for QRPA and IBM would produce a narrower 
distribution when qa is incorporated with the matrix elements. 

It could be anticipated that some correlation between the energy release in a 
decay and the size of the matrix element might arise simply from differences in 
the nuclear structure between parent and progeny. In the shell model, a large Q- 
value and high Z both tend to accentuate differences in the neutron and proton 
valence shells involved in the transition. However, because of the spin-orbit force, 
a monotonic trend is not in fact expected. Specifically, ^"^Mo in the middle of the 
suite of isotopes is one of the few cases (perhaps the only case) where the spin-orbit 
partners are open in the ground states: the ^§7/2 valence neutrons can decay into 
TTgg/2 proton states. A much larger space is needed for a good description, naturally, 
and there is unfortunately no large-basis shell model calculation for this nucleus. 
Nevertheless, the QRPA, GCM, IBM, and PHFB models tend to give somewhat 
larger values for ^"^Mo, possibly reflecting this favorable situation. In the IBM, the 
matrix elements display a parabolic dependence that peaks between closed shells and 
have a general dependence as A~'^^^^ A comprehensive understanding, across the 
theoretical spectrum, of the relationship captured in Eq. [4] between the phase space 
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and the matrix element could provide useful additional insight into the predictive 
accuracy of the calculations of the nuclear matrix elements. 

3. Conclusions 

Large-scale experiments searching for neutrinoless double beta decay are in the plan- 
ning and prototype stages. The importance of the objective has led to a productive, 
intensive, and world-wide theoretical effort to make the best possible predictions 
for the nuclear matrix elements. That effort continues apace. 

We have organized currently available theoretical input on nuclear matrix ele- 
ments, phase-space factors, and gA in a way that is intended to be helpful to ex- 
perimental design teams and to provide feedback to theoretical teams. The results 
are somewhat surprising. First, we find that there is little evidence on theoretical 
grounds to favor or disfavor the choice of one isotope over another. All eleven for 
which a significant body of theoretical calculations exists seem to have about the 
same sensitivity to Qvfjj3 decay per unit mass. Second, we note a striking inverse 
correlation between the phase-space available and the size of the nuclear matrix 
element. There are essentially no exceptions, not even ^®Ca, provided the Q-value 
is > 2 MeV. Third, barring a large theoretical bias common to all existing theo- 
retical methods, experiments at the tonne scale with negligible background can be 
expected to be sensitive to Majorana masses in the 10-100 meV range. 

The signal reported by Klapdor-Kleingrothaus and Krivosheina^^^ for QvjSjS in 
'^^Ge is 246(40) decays per tonne per year (1 tonne = 1 Mg = 1000 kg). Upper limits 
have recently been reported for Oz^/3/3 in ^'^^Xe by the EXO Collaboration^^ and the 
KamLAND-Zen Collaboration.'^^ The results are summarized in Table H] The most 



Table 2. Comparison of recent results on the rate for 
Ov/3/3 decay. 



Isotope 


Specific Activity 


Reference 




Mg-l y-i 




76Ge 


246 ± 40 


[16] 


"6Xe 


< 192 (90% CL) 


en 


i36Xe 


< 162 (90% CL) 


m 


Combined: 






i36Xe 


< 90 (90% CL) 





probable value for the difference between two samples drawn at random from a nor- 
mal distribution is one standard deviation. Within the framework of the empirical 
scaling rule Eq. [4] the difference in specific activity expected from the distribution 
in matrix elements is about a factor of 3. (The uncertainty in the renormalization of 
gA is omitted for this purpose on the premise that the effective value of qa is likely 
to have much in common from isotope to isotope.) As it happens, the experimental 
rate limit from ^'^^Xc differs by about a factor of 3 from the positive result reported 
for ^^Ge, a not improbable outcome in this analysis. Although the -'^^^Xe limit is a 
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90% CL upper limit rather than a positive observation, the extent of disagreement 
is not substantial. The difficulty in drawing conclusions from results in two differ- 
ent isotopes comes from the width of the distribution of specific activities (Fig. [s]) 
arising from theoretical variance, which is nearly normal in the logarithm. A more 
detailed statistical analysis without using the universality relation outlined here has 
been presented recently by Bergstromj^i^l who reaches the conclusion that the Ge 
and Xe results are probably incompatible. However, the widths of the individual 
nuclear matrix element distributions assumed are much narrower. 

As experiments searching for 0i//3f3 advance toward very large scales, the ab- 
solute specific activity appears to be largely independent of the choice of isotope. 
Instead, the decision on which isotopes offer the best opportunities for scale-up 
may be based on more technical criteria: cost, redundant identification of candidate 
Oz//?/3 events, background, resolution, and similar issues. We conclude first that a 
theoretical understanding of the correlation between the nuclear matrix element and 
the effective value of the axial- vector coupling constant gA offers perhaps the most 
immediate path toward reducing the width of the distribution shown in Fig. [5] Sec- 
ondly, theoretical understanding of the inverse correlation between the phase space 
and the square of the nuclear matrix element will add confidence in the predictive 
accuracy of the theory of neutrinoless double beta decay. 
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